Magnetic Anisotropy Variations and Non-Equilibrium Tunneling in a Cobalt 

Nanoparticle 
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We present detailed measurements of the discrete electron-tunneling level spectrum within 
nanometer-scale cobalt particles as a function of magnetic field and gate voltage, in this way prob- 
ing individual quantum many-body eigenstates inside ferromagnetic samples. Variations among 
the observed levels indicate that different quantum states within one particle are subject to difi^er- 
ent magnetic anisotropy energies. Gate-voltage studies demonstrate that the low-energy tunneling 
spectrum is affected dramatically by the presence of non-equilibrium spin excitations. 

PACS numbers; 73.22.-f, 73.23.Hk, 75.75.+a 



The electronic states within ferromagnetic devices are 
understood surprisingly poorly. For example, several 
types of experiments have been done to characterize spin 
polarization near the Fermi level , but different tech- 
niques give different results. Part of the difficulty is 
that these experiments average over large numbers of 
states, and different types of experiments effectively take 
differently weighted averages. Recently, we have devel- 
oped a tunneling technique that can resolve the individ- 
ual states within nm-scale samples of ferromagnets 
This method has helped motivate new descriptions of 
ferromagnetism that go beyond mean-field Stoner models 
[@,|| . However, the first experiments left many open ques- 
tions, particularly concerning the proper description of 
anisotropy energies, and whether the tunneling spectrum 
refiects the true electronic density of states or whether it 
is modified by non-equilibrium effects. To answer these 
questions, we have introduced a gate electrode to our 
devices and have also developed lower-noise tunnel bar- 
riers, thereby allowing new types of measurements and 
greatly improving the data quality. We present evidence 
that the quantum states in a ferromagnetic nanoparticle 
are not all described by the same anisotropy-energy func- 
tion that governs the ground state |^ , but the anisotropy 
varies from state to state. We show that non-equilibrium 
processes induced by tunneling affect the measured spec- 
trum, thereby explaining a larger-than anticipated den- 
sity of resonances. The gate electrode also allows a com- 
parison with recent theories that minority-electron 
tunneling should dominate in a nanoparticle. 

We will show data from both gated and non-gated de- 
vices. In either case, the samples contain a Co nanopar- 
ticle, separated from aluminum electrodes by aluminum 
oxide tunnel barriers, inside a tunnel junction with a 
nano-scale area small enough to allow individual parti- 
cles to be contacted. (See device schematic. Fig. 1(a).) 
The devices without a gate electrode are fabricated 
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FIG. I. (a) Cross-sectional device schematic, (b) and 
(c) Colorscale plots of dl/dV for tunneling resonances in a 
cobalt nanoparticle. The field is varied from positive to neg- 
ative values, and dl /dV vs. V is measured at each step. The 
maximum conductance is 3 x 1Q~''^Q,~^ . 

using a procedure described previously , with the new 
innovation that by using 1 nm of deposited aluminum 
oxide to form the tunnel barrier on the cobalt particle 
we produce junctions with reduced charge noise. The 
nanoparticles are made by depositing 0.5 nm of Co at 
room temperature, which makes particles in the range 
of 1—4 nm diameter The gated device is made by 
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forming a hole about 10 nm in diameter in a suspended 
silicon-nitride membrane, depositing 18.5 nm of Al to 
make the gate electrode on the lower side of the device 
as shown in Fig. 1(a), and then isolating this gate by an- 
odizing in an oxygen plasma to 3.5 V bias and depositing 
8.5 nm of SiOx. The rest of the fabrication proceeds in 
the same way as for the non-gated devices. 

When cooled to dilution refrigerator temperatures, the 
tunnehng conductance dl/dV as a function of source- 
drain voltage V consists of individual peaks, associated 
with transitions between discrete electronic states in the 
nanoparticle Q. The T/-spacings of resonances can be 
converted to energy, = e/\VC2/ {Ci -\- C2), in this 
way correcting for capacitive division across the two junc- 
tions. The capacitance ratio can be determined by com- 
paring peak positions at positive and negative V |p^ , |rT| . 
In Fig. l(b,c), we plot the energies of tunneling reso- 
nances for a Co nanoparticle in a non-gated device as a 
function of magnetic field, H . As H is swept from posi- 
tive values toward zero, the levels first undergo significant 
continuous shifts. The discontinuity near = is an arti- 
fact of the Al leads going superconducting and then being 
driven normal by a negative field. This causes the reso- 
nance energies to jump by A, the superconducting gap 
, but the energies of the states within the nanopar- 
ticle evolve continuously. Near fioHsw = —0.120 T, all 
the levels exhibit another large discontinuity, which can 
be identified with magnetic switching of the nanoparti- 
cle li H is swept from negative to positive values, 
the field value for this transition is hysteretic, and the 
conductance spectrum is identical to that shown, but re- 
flected about H=0. At larger values of \H\ (Fig. 1(c)), 
all the levels move to increasing energy, with slopes that 
correspond to effective (7-factors ranging from 0.06 to 1.1. 
The fact that the ^-factors are reduced below 2 indicates 
that the resonances are not purely spin-up or spin-down, 
but are mixed by spin-orbit scattering [ |l2[ . 

These data confirm in more detail several observations 
made previously |6|: the mean energy spacing is about 
0.2 meV, smaller than estimated for independent elec- 
trons within a particle 1—4 nm in diameter; the energies 
are coupled to the orientation of the particle's magnetic 
moment; and resonances with a given high-field slope can 
jump either up or down at Hsw In addition, these data 
reveal several important new features. (1) More reso- 
nances are resolved than in the previous work, and it is 
clear that every level has a different H dependence. (2) 
When the resonances approach close to each other as a 
function of H, they cross, with no evidence of avoided- 
crossing behavior. (3) The levels can exhibit complicated 
non-monotonic variations as a function of H , for instance 
the levels marked with dots in Fig. 1(b). 

The new observations are incompatible with the sim- 
plest model proposed previously for states in magnetic 
nanoparticles Q, which assumed that all the states in 
a particle could be described by a single anisotropy- 
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FIG. 2. Tunneling energy AEfi calculated using the semi- 
classical model discussed in the text, with H at 45° from the 
easy axis, for various values of the anisotropy-energy differ- 
ence Sk — fc]v+i — fcjv. In (a) the total spin is increasing; in 
(b) it is decreasing. The curves are offset for clarity. 

energy function. This assumption has the consequence 
that all transitions in which the spin component (Sz) in- 
creases should exhibit the same H dependence, the (Sz)- 
decreasing transitions should show a second behavior, 
and neither dependence should exhibit the complicated 
non-monotonicities below 1 Tesla evident in Fig. 1(b). 
(The results of the older model are shown in the Sk — 
curves in Fig. 2.) The model in |^ assumed a uniaxial 
anisotropy, but we have confirmed that the same behav- 
iors hold for more complicated forms which include terms 
of 2nd and 4th-order in the spin-components |^ , as long 
as one assumes that the same anisotropy function applies 
to all electronic states. 

The fact that all the tunneling resonances undergo dif- 
ferent energy variations in the \ow-H range where the 
magnetic moment is being reoriented indicates directly 
that all the electronic states of the particle cannot be 
described by the same anisotropy-energy function. We 
have explored whether such variations may also affect the 
form for the 77-dependence of the energies as described in 
point (3) above, and we find that they provide a natural 
explanation for the complicated non-monotonic behavior 
as a function of H. We start by considering one sin- 
gle resonance associated with a transition between two 
states with N and iV -t- 1 electrons. We extend the N- 
electron Hamiltonian stated in |^ in the simplest way to 
incorporate variations in anisotropy energy: 

1-1 = -gcsfJ'B^J-oH ■ S ~kN{S -71)^/30. (1) 

Here hS is the total spin with ground-state magnitude 
hSa for N electrons, n is a unit vector in an easy-axis 
direction, and the (uniaxial) anisotropy energy prefactor 
kjsi is now allowed to vary between the N- and (iV-|-l)- 
electron states. For simplicity, we assume that the easy 
axis is the same for all states. We have solved for the 
ground state energies for N and N + 1 electrons semi- 
classically as a function of H by finding the spin ori- 
entation that gives the local minimum in Eq. (1), as- 
suming that So does not vary with H M, and then we 
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calculate the form of the tunneling transition energies as 
E{N+1, H) - E{N, H). The results are shown in Fig. 2. 
Assuming Sq ^ 1000 (appropriate to a 4 nm Co parti- 
cle), and an average value of (fcAr) « 0.01 meV in accord 
with the switching field [pp^, fluctuations in /cat of order 
1-3% are sufficient to explain both the size and form of 
the non-monotonicities. Subsequent to this initial semi- 
classical analysis, similar conclusions were also reached 
in a more rigorous quantum-mechanical picture ^ . 

Although fluctuations in the properties of eigenstates 
are not often considered in the context of ferromagnets, 
they are not surprising. In non-magnetic particles, the 
g-factors for Zeeman splitting fluctuate and the 

statistics of these fluctuations have been investigated in 
random-matrix treatments of the spin-orbit interaction 
p^ , p^ . Anisotropy in magnetic particles also arises from 
spin-orbit interactions. An additional consequence of 
anisotropy-energy fluctuations should be that the value 
of Hsw will vary by 1-3% depending on the occupation of 
excited electronic states. This has yet to be investigated. 

We now turn to data from a gated tunneling device, 
with the primary motivation being to resolve the ques- 
tion 1^ of whether the larger-than-anticipated density 
of tunneling resonances observed at low energy (noted 
above) is due to inelastic emission of spin-waves during 
tunneling |16| or due to non-equilibrium effects |p|,p7t. 
The idea of the non-equilibrium mechanism is that the 
energy of tunneling transitions can be described as the 
energy differences of states in the nanoparticle with N 
and TV ± 1 electrons: l^Ef^ = Ef^^ - E^ . If under 
conditions of current flow, E^ can assume an ensemble 
of different values, for instance due to low-energy spin 
excitations or electron-hole excitations within the parti- 
cle, then the number of observed tunneling resonances 
can increase above the number that originate only from 
a single equilibrium ground state. We can test this with 
a gate voltage because the energy of a tunneling transi- 
tion can be tuned from high values down close to zero 
where tunneling can be initiated by small V . For suf- 
ficiently small V , the tunneling electrons may have in- 
sufficient energy to excite non-equilibrium states. There- 
fore, a test of whether tunneling resonances are associ- 
ated with non-equilibrium initial states is whether some 
transitions disappear when they are tuned to small V . 
This disappearance is exactly what is observed for the 
lowest-energy transitions at 7?=0 (Fig. 3(a)), when the 
electrodes are superconducting. 

Related non-equilibrium effects have been observed 
in non-magnetic particles, but the consequences are 
much less dramatic. In Al particles, the energy shifts 
due to non-equilibrium are small, resulting only in a 
fine-structure about the energy of equilibrium tran- 
sitions |jl^,|l8|. We have not observed well-resolved 
transitions in non-magnetic particles to disappear com- 
pletely as Vg is varied, whereas at least the first 5 
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FIG. 3. Colorscale conductance plots of a gated Co 
nanoparticle as a function of Vg and V . (a) Superconducting 
leads (-^=0). (b) Normal-state leads (0.07 Tesla). The dashed 
lines in (a) indicate the expected evolution of the threshold 
peaks if the resonances did not involve non-equilibrium pro- 
cesses. The voltage resolution is better at H=0 because of the 
superconducting singularity in the electrode density of states. 
The electron T f« 90 mK. 

lowest-energy states in the Co particle lose conductance. 
Non-equilibrium effects therefore appear to be much 
stronger in Co, perhaps due to larger fluctuations in 
electron-electron interactions and/or a large multiplicity 
of low-energy collective spin excitations (in addition to 
the particle-hole excitations that were considered previ- 
ously [p^ ). Since the relaxation rate of non-equilibrium 
excitations must be slower than the tunneling rate for 
the spectrum to be affected, the relaxation rate is of the 
order of or slower than ^ I MHz. 

The presence of level crossings (instead of avoided 
crossings), noted in Fig. 1 provides independent new ev- 
idence supporting the non-equilibrium scenario. In non- 
magnetic particles, when spin-orbit scattering reduces 
the large- g- factors to less than 1.7, tunneling transi- 
tions originating from the same initial state exhibit clear 
level repulsion p2| . In contrast, under non-equilibrium 
conditions, tunneling resonances occurring at similar val- 
ues of V can result from different pairs of eigenstates 
{E^ , E^^^), none of which are nearly degenerate, so an 
avoided crossing would not be expected. 

Despite these two lines of evidence for the importance 
of non-equilibrium transitions, our observations are not 
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FIG. 4. High-field conductance spectra for the sample of 
Fig. 3 at (a) Vg = 20 mV and (b) Vg = 125 mV. 

in full agreement with the simplest phenomenological 
scenario that includes non-equilibrium effects Ref. 

proposes that non-equilibrium spin-accumulation may 
occur by a sequence of transitions up a ladder of energy 
states having different total spin S*, with the sequence ter- 
minating when the next step up the ladder requires more 
energy than that provided by the source-drain bias, eV . 
Within this scenario, a fraction of both the low-energy 
and high-energy tunneling resonances should disappear 
as Vg tunes the transitions to lower voltages. In Fig. 3(a), 
only the lowest-energy transitions experience a loss of 
tunneling amplitude. We suspect that a more compli- 
cated set of states than was used in may need to be 
considered in characterizing non-equilibrium processes on 
our magnetic nanoparticle. For instance, single-electron 
quasiparticle excitations and also ladders of states with 
the same value of total spin S but different values of Sz 
may be involved. 

Another unanticipated observation is that the disap- 
pearance of tunneling resonances as a function of Vg oc- 
curs only when the electrodes are superconducting, and 
not when they are normal (Fig. 3(b)). Scans at 0.07, 1, 
3, 5, and 8.5 T are all qualitatively similar. Perhaps, be- 
cause of the singularity in the density of states for super- 
conducting electrodes, the ensemble of non-equilibrium 
states may be populated differently by tunneling from 
normal and superconducting electrodes, with less non- 
equilibrium at low V in the superconducting case. How- 
ever, this process has not yet been modeled in detail. 

One more simple benefit of having a gate electrode is 
that it allows a determination of whether the tunneling 
resonances correspond predominantly to the motion of 
majority or minority electrons. Consider the data shown 
in Figures 3 and 4. At = 20 mV the low-energy tun- 
neling thresholds correspond to the addition of an elec- 
tron to the nanoparticle {N ^ N +\), and at large mag- 
netic fields (Fig. 4(a)) the tunneling energies increase, 
meaning that the tunneling states have lower {Sz)- It 
follows that the transitions correspond predominantly to 
minority-electron tunneling. Similarly, at Vg = 125 mV 
(Fig. 4(b)), the tunnehng thresholds are A^-|-1^A^ pro- 



cesses in which the spin increases, so a minority electron 
is tunneling off the nanoparticle. These results from one 
sample are in agreement with predictions [0J|] that the 
total spin should almost always decrease when electrons 
are added, due to both a large minority density of states 
and strong exchange interactions. Of course, because 
of spin-orbit coupling, the transitions are not purely of 
minority-electron character. 

In summary, measurements of the TJ-dependence of 
the tunneling-resonance energies in cobalt nanoparticles 
exhibit level-to-level differences and non-monotonic vari- 
ations that can both be explained by variations of 1-3% 
in magnetic anisotropy energy between different states. 
The observation of level crossings and the disappearance 
of resonances as a function of Vg indicate that the major- 
ity of low-energy resonances are associated with tunnel- 
ing transitions from non-equilibrium initial states. The 
strength of non-equilibrium effects appears to depend on 
whether the electrodes are normal or superconducting. 
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